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(57) A signal component (e.g. , a period signal com- 
ponent such as a sinusoid) is added (21 0) to an original 
color image to form a modified color image, with the goal 
of embedding the signal in the image so that it is imper- 
ceptible to a human viewer. A comparing operation 
(240) uses a model of human perception to measure the 
perceptual difference between the original and modified 
images, identifying local areas of the modified image 
where the signal difference exceeds a threshold, indi- 
cating that the signal is perceptible to a human viewer. 
Using the perceptual difference measurement data; the 
signal is attenuated in the identified local areas that in- 
dicate a perceptually unacceptable difference, and this 
modified signal component is then added (210) to the 
original color image in a next iteration. Perceptual dif- 
ference measurement and signal attenuation are iterat- 
ed until the comparison operation determines that the 
signal difference is perceptually acceptable. The tech- 
nique takes advantage of the fact that, although the spa- 
tial frequencies of the embedded signals are well within 
the range of spatial frequencies to which humans are 
normally quite sensitive in the luminance (black-white) 
vision channel, this sensitivity does not extend to the 
color vision bands. In an illustrated embodiment, a set 
of sinusoidal signals that form a grid are added to the 
color image; location (i.e., decoding) of the sinusoids, 
which does not require the original color image, allows 
computing a geometric mapping from an image with the 
embedded signals to the original image. 
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Description 

[0001] Encoding information in image form to permit its subsequent electronic decoding is a well-known information 
processing technique. For example, bar codes explicitly carry encoded information in black and white image form, and 
5 are typically used in applications where the obvious and perceptible presence of the encoded information is intended 
and is not a disadvantage. 

[0002] Data glyph technology is a category of embedded encoded information that is particularly advantageous for 
use in image applications that require the embedded data to be robust for decoding purposes yet inconspicuous, or 
even surreptitious, in the resulting image. Data glyph technology encodes digital information in the form of binary 1's 
10 and 0's which are then rendered in the form of distinguishable shaped marks such as very small linear marks. U.S. 
Patent 5 : 315 : 098 : entitled "Methods and Means for Embedding Machine Readable Digital Data in Halftone Images" 
discloses techniques for encoding digital data in the angular orientation of circularly asymmetric halftone dot patterns 
that are written into the halftone cells of digital halftone images. 

[0003] Another type of image or document marking is known as digital watermarking. A successful digital watermark- 
15 ing technique simultaneously achieves two purposes: first, the technique must produce an embedded signal that is 
imperceptible to a human viewer so as not to diminish the commercial quality and value of the image being watermarked. 
At the same time : the embedded signal must be resistant to tampering; removal of the embedded signal defeats the 
identification purpose of watermarking, and so a successful watermarking technique is typically designed so that at- 
tempts to remove the embedded signal cause degradation of the image sufficient to render it commercially less valuable 
20 or worthless. 

[0004] Brassil et al., in "Electronic Marking and Identification Techniques to Discourage Document Copying" in IEEE 
Journal on Selected Areas in Communications, Vol. 1 2, No. 8, October 1 995, pp. 1 495 - 1 504, disclose three techniques 
for embedding a unique codeword in a text document image that enables identification of the sanctioned recipient of 
the document while being largely indiscernible to document readers, for the purpose of discouraging unauthorized text 
2S document distribution. The image coding schemes were designed to result in a substantial loss of document presen- 
tation quality il successfully removed. The techniques disclosed include line shift coding, word shift coding and feature 
coding. Use of these techniques in the resulting image is typically not noticeable to a viewer of the image, and text in 
the image is not substantively altered. 

[0005] Other techniques are described in WO-A-95/14289, US-A-5278400 and Cox, Kilian, Leighton and Shamoon, 

30 in NEC Research Institute Technical Report No 95 10 entitled "Secure Spread Spectrum Watermarking for Multimedia". 
[0006] Many of the existing techniques for embedding information in images appear to operate in the black and white 
image domain, and so do not explicitly address how to embed a signal in a color image that is imperceptible to a human 
viewer and that does not distort the quality of the image. Digital watermarking techniques, even those that may apply 
to color images, are typically designed to be irreversible; they produce a tamper-proof embedded signal which cannot 

35 be removed without distorting the information in the image; the watermarked image must remain watermarked for all 
subsequent uses. Moreover, the detection of an embedded identification signal in a watermarked image typically re- 
quires the use of the original image, which is typically maintained in a secure location for such future use as needed. 
While these characteristics of digital watermarking are useful features for image authentication and identification pur- 
poses, they may be limitations for other purposes. 

40 [0007] In accordance with a first aspect of the present invention, we provide a method for operating a machine to 
produce a version of an original color image data structure having perceptually acceptable signal differences embedded 
therein: the machine including a processor and a memory device for storing data; the data stored in the memory device 
including instruction data the processor executes to operate the machine: the processor being connected to the memory 
device for accessing and executing the instruction data stored therein; the method comprising: operating the processor 

45 to obtain an original colorimage data structure defining an original color image including a plurality of original image 
locations each indicating an original color; operating the processor to obtain a signal; operating the processor to add 
the signal to the original color image data structure to produce an interim color image data structure: the interim color 
image data structure being a version of at least a portion of the original color image data structure and including a 
plurality of interim image locations each indicating an interim color: the signal being added to at least one of the plurality 

so of interim image locations: an original image location in the original color image data structure and an interim image 
location in the interim color data structure indicating the same image location in each image and being referred to 
collectively as respectively paired image locations, each one of the respectively paired image locations indicating an 
image color; operating the processor to perform a comparison operation to produce a quantitative perceptual difference 
measurement at each of the respectively paired image locations; the comparison operation using a model of human 

55 perception to quantitatively represent the perceptual differences as perceived by a human viewer of the original and 
interim color images; operating the processor to attenuate the signal to produce a non-zero modified signal when the 
difference measurements between the respectively paired image locations exceed a threshold indicating an acceptable 
perceptual difference between the original and interim color images; and operating the processor to perform an iteration 
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operation iterating the adding, comparison and attenuating operations using the modified signal until the perceptual 
difference measurements between the respectively paired image locations do not exceed the threshold; the iteration 
operation producing a modified color image data structure indicating a perceptually acceptable signal difference be- 
tween the original and modified color images. 

5 [0008] In accordance with a second aspect of the present invention, we provide a method for operating a machine 
to add humanly imperceptible information to an original color image data structure to produce a modified color image; 
the machine including a processor and a memory device for storing data; the data stored in the memory device including 
instruction data the processor executes to operate the machine: the processor being connected to the memory device 
for accessing and executing the instruction data stored therein; the method comprising: operating the processor to 

to obtain an original color image data structure defining an original color image including a plurality of original image 
locations each indicating an original color; operating the processor to obtain periodic signals indicating information; 
operating the processor to add the periodic signals indicating information to the original color image by modifying a 
plurality of the original image colors in original image locations to produce an interim color image data structure; the 
interim color image data structure being a version of at least a portion of the original color image data structure and 

*5 including a plurality of interim image locations each indicating an interim color: an original image location in the original 
color image data structure and an interim image location in the interim color data structure indicating the same image 
location in each image and being referred to collectively as respectively paired image locations; each one of the re- 
spectively paired image locations indicating an image color; operating the processor to perform a comparison operation 
to produce a quantitative perceptual difference measurement at each of the respectively paired image locations; the 

20 comparison operation using a model of human perception to quantitatively represent the perceptual differences as 
perceived by a human viewer of the original and interim color images: operating the processor to attenuate the periodic 
signal to produce a non-zero modified periodic signal when the difference measurements between the respectively 
paired image locations exceed a threshold indicating an acceptable perceptual difference between the original and 
interim color images; and operating the processor to perform an iteration operation iterating the adding, comparison 

25 and attenuating operations using the modified signal until the perceptual difference measurements between the re- 
spectively paired image locations do not exceed the threshold; the iteration operation producing a modified color image 
data structure including the humanly imperceptible information and indicating a perceptually acceptable signal differ- 
ence between the original and modified color images. 

[0009] In accordance with a third aspect of the present invention, we provide a method for operating a machine to 

30 automatically transform geometric properties of an acquired image version of an original image to match geometric 
properties of the original image; the machine including a signal source connected for receiving signals, a processor 
and a memory device for storing data; the data stored in the memory device including instruction data the processor 
executes to operate the machine; the processor being connected to the memory device for accessing and executing 
the instruction data stored therein; the method comprising: operating the processor to obtain from the signal source 

35 an acquired image data structure defining an acquired image: the acquired image being a version of an original image 
and additionally having embedded signals therein not included in the original image; the embedded signals having 
predetermined geometric relationships with respect to each other; operating the processor to obtain geometric con- 
straint data indicating expected geometric relationships about the embedded signals in the acquired image: operating 
the processor to determine the predetermined geometric relationships of the embedded signals in the acquired image; 

40 operating the processor to compute geometric differences between the acquired image and the original image using 
the predetermined geometric relationships of the embedded signals in the acquired image and using the geometric 
constraint data indicating the expected geometric relationships of the embedded signals; operating the processor to 
transform geometric properties of an acquired image version of an original image to match geometric properties of the 
original image using the geometric differences computed between the acquired and original images. 

45 [0010] In accordance with a fourth aspect of the present invention, we provide a method for operating a processor- 
controlled machine to decode periodic signals embedded in an image; the machine including a processor and a memory 
device for storing data; the data stored in the memory device including instruction data the processor executes to 
operate the machine; the processor being connected to the memory device for accessing and executing the instruction 
data stored therein, the method comprising: operating the processor to obtain an image data structure defining an 

50 image: the image including a set of n periodic signals embedded therein; the set of n periodic signals having geometric 
relationships with respect to each other; operating the processor to obtain geometric constraint data indicating expected 
geometric relationships about the set of n periodic signals embedded in the image; operating the processor to identify 
a plurality of local peak power concentrations in a two-dimensional spatial frequency power spectrum of the image: 
operating the processor to eliminate a subset of the local peak power concentrations that do not satisfy the geometric 

55 constraint data indicating expected geometric relationships about the set of n periodic signals embedded in the image; 
a remaining subset of the local peak power concentrations being determined to be candidate local peak power con- 
centrations indicating spatial frequencies and directions of the set of n periodic signals; and operating the processor 
to compute a best linear mapping between sets of n candidate local peak power concentrations in the image and a set 
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of n periodic signals indicated by the geometric constraint data; the best linear mapping identifying the spatial frequen- 
cies and directions of the set of n periodic signals embedded in the image. 

[0011] The steadily rising use of color images in all types of commercial and aesthetic applications suggests that 
many of the techniques that have been developed for embedding information in black and white images need to be 

5 extended to the color domain. The present invention is premised on the observation that modifying a color image by 
adding signals that do not unacceptably distort or degrade the image to a human viewer presents a human perception 
problem that is different from that associated with black and white images. The mere extension of existing techniques 
in the black and white image domain using empirical or heuristic approaches to determine whether an added signal is 
humanly perceptible is inadequate to embed signals in color images that do not unacceptably distort or degrade the 

to image. Moreover, image artifacts that appear in a color image as a result of adding a signal using an inadequate method 
are likely to be attributed to the hardware device that produced the image; if the quality of the color image is aesthetically 
unacceptable, the hardware device or the application that produced these artifacts will simply not be used. 
[001 2] The present invention is premised on the discovery that adding signals to a color image that do not unaccept- 
ably distort or degrade the image can only be predictably and reliably accomplished by using a sophisticated model 

is of human perception that is able to quantitatively determine the magnitude of the perceptual differences between an 
original image and a version of the image with the embedded signals. The technique of the present invention uses a 
quantitative model of human perception to attenuate the power (amplitude)of the added signal in local regions of the 
color image where the model indicates that the perceptual difference between an original color and the modified color 
produced with the added signal is too high, toward the goal of producing a version of an original color image having 

20 an added signal that is substantially imperceptible to human viewers of the image. The quantitative model of human 
perception controls the perceptibility of the embedded signal by ensuring that it is below a perceptual threshold. 
[0013] The technique of an illustrated embodiment of the present invention adds a pattern of periodic signals to a 
color image, and in particular, adds a pattern of amplitude-modulated sinusoidal signals to the color image. A comparing 
operation automatically determines local areas in the modified color image where the amplitude of the embedded 

2S information is too high : and is thus perceptible. The perceptible signals are then iteratively attenuated in those identified 
local areas of the modified image. In one aspect of the illustrated embodiment, the sinusoidal signals are added to the 
yellow-blue opponent-color band of the color image, at spatial frequencies where most color images have relatively 
tittle power, and where humans have the least sensitivity. The added amplitude-modulated sinusoidal signals can be 
decoded (located in a modified image) because they form a specific pattern of peak power concentrations in the fre- 

30 quency domain. The manner of adding the signals results in particular geometric relationships occurring between the 
spatial frequencies that are unlikely to occur by chance in natural or computer-generated synthetic images. The em- 
bedding technique takes advantage of human perception: the spatial frequencies of the embedded sinusoidal signals 
are well within the range of frequencies to which humans are normally quite sensitive in the luminance (black-white) 
vision band, but this sensitivity does not extend to the color vision bands. Thus, while sinusoidal signals at relative high 

35 spatial frequencies are added to the modified color image in the embodiment of the invention described below, signals 
can be added at tower spatial frequencies if a particular application, a specific decoding domain, or a specific device 
requires it. 

[0014] The embedded sinusoidal signal may be useful in a variety of applications. A set of sinusoids form a grid 
when embedded in an original image. After a decoding operation has extracted the spatial frequencies of the embedded 

40 sinusoids, a mapping is computed between the acquired (e.g., scanned) modified image and the sinusoidal grid image 
that was embedded in the original image. This mapping then allows geometric transformations to be performed on the 
acquired image to match the geometric properties of the acquired image to the geometric properties of the original 
image. For example, the acquired image can be automatically aligned and sealed to that of the original image, if the 
acquired image has become rotated, skewed, reduced or enlarged during previous manipulations. Thus, images con- 

45 taining the added sinusoidal signals do not have to be perfectly aligned when scanned, nor manually sealed or rotated 
after scanning. In addition, the sinusoidal grid can function as a reference point with respect to which other embedded 
information can be located with precision. Moreover, a combination of embedded sinusoidal signals can be chosen 
such that the total embedded signal has a period greater than or equal to the image field so that, during a decoding 
operation, each position in the acquired image is uniquely associated with a unique position in the embedded signal 

so image. 

[0015] A significant advantage of the technique of the signal encoding process is that the process of locating the 
sinusoidal signal image that is embedded in an acquired image (i.e., the decoding process) does not require use of 
the original, unmodified image. This characteristic provides anyone with the ability to use the embedded information. 
[0016] Therefore, in accordance with one aspect of the present invention, a method is provided for operating a proc- 
55 essor-controlled machine to produce a version of an original image data structure having perceptually acceptable signal 
differences embedded therein. The machine includes a signal source for receiving data; memory for storing data; and 
a processor connected for accessing and executing instruction data stored in the memory for operating the machine; 
the processor is further connected for receiving data from the signal source; and connected for storing data in the 
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memory. The method comprises operating the processor to obtain an original image data structure defining an original 
image that includes a plurality ot original image locations each indicating an original color, and to obtain a signal for 
embedding in the original image The method further comprises operating the processor to add the signal to at least 
one of the plurality of image locations of the original image data structure to produce an interim color image data 
s structure. The interim color image data structure is a version of at least a portion of the original image data structure 
and includes a plurality of interim image locations each indicating an interim color. An original image location in the 
original image data structure and an interim image location in the interim color data structure indicate the same image 
location in each image and are referred to collectively as respectively paired image locations. Each one of the respec- 
tively paired image locations indicate an image color. A comparison operation is then performed to produce a quanti- 
se tative perceptual difference measurement at each of the respectively paired image locations. The comparison operation 
uses a model of human perception to quantitatively represent perceptible image differences as perceived by a human 
viewer of the original and interim images. When the difference measurements between the respectively paired image 
locations exceed a threshold that indicate an acceptable perceptual difference between the original and interim color 
images, the signal is attenuated to produce a non-zero modified signal. Then an iteration operation is performed that 
75 iterates the adding, comparison and attenuating operations using the modified signal until the perceptual difference 
measurements between the respectively paired image locations do not exceed the threshold. The iteration operation 
produces a modified color image data structure that indicates a signal difference between the original and modified 
color images that is a perceptually acceptable signal difference. 

[0017] In accordance with another aspect of the invention, a method for operating a processor -control led machine 
20 js provided for decoding a set of n periodic signals embedded in an image. The method comprises obtaining an image 
data structure defining an image including a set of n periodic signals embedded therein. The set of n periodic signals 
have geometric relationships with respect to each other. The method further includes obtaining geometric constraint 
data indicating expected geometric relationships about the set of n periodic signals. Then, a plurality of local peak 
power concentrations in the two-dimensional spatial frequency power spectrum of the image are identified; a subset 
25 of these local peak power concentrations that do not satisfy the geometric constraint data indicating the expected 
geometric relationships about the set of n periodic signals are eliminated from consideration. A remaining subset of 
the local peak power concentrations are determined to be candidate local peak power concentrations indicating spatial 
frequencies and directions of the set of n periodic signals. Then, a best linear mapping is computed between sets of 
n candidate local peak power concentrations in the image and a set of n periodic signals indicated by the geometric 
30 constraint data The best linear mapping identifies the spatial frequencies and directions of the set of n periodic signals 
embedded in the image. 

[0018] Some examples of methods and apparatus according to the present invention will now be described with 
reference to the accompanying drawings, in which:- 

35 Figure 1 is a flowchart illustrating the general operation of the invention for adding a signal component to a color 

image according to the invention; 

Figure 2 schematically illustrates a representative original color image that is suitable for input to the present 
invention; 

Figure 3 schematically illustrates a representative signal component image that is suitable as the signal component 
40 to be added to the original color image of Figure 2 according to the present invention; 

Figure 4 schematically illustrates a modified version of the original color image of Figure 2 with the added signal 
component image of Figure 3; 

Figure 5 schematically illustrates a perceptual difference image showing the perceptual differences between the 
original image of Figure 2 and the modified color image o1 Figure 4: 
45 Figure 6 schematically illustrates the signal component image of Figure 3 modified with an attenuated signal com- 

ponent in certain image regions: 

Figure 7 is a block diagram schematically illustrating the general operation of an illustrated embodiment of the 
present invention; 

Figure 8 is a flowchart illustrating the general operation of the S-CIELAB color image comparison operation of the 
50 illustrated embodiment of Figure 7; 

Figure 9 is a flowchart illustrating the technique for producing the attenuation image according to the illustrated 
embodiment of Figure 7; 

Figure 10 schematically illustrates an image of a set of sinusoidal signals that are suitable for embedding in a color 
image according to an illustrated embodiment of the present invention; 
55 Figure 11 schematically illustrates a cross section of the image of Figure 10 showing a pair of the sinusoids in that 

image; 

Figure 1 2 schematically illustrates a second image of a set of sinusoidal signals that are also suitable for embedding 
in a color image according to the present invention; 
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Figure 1 3 is a flowchart that illustrates the operations of a class of applications which uses an image having signals 
embedded according to the technique of the present invention; 

Figure 14 is a flowchart of a decoding operation for locating the embedded signal in an acquired image according 
to an illustrated embodiment; 

s Figure 15 is an amplitude spatial frequency spectrum image of the yellow-blue band of a modified color image 

having the signals of Figure 10 embedded therein, and showing certain geometric relationships between the peak 
power concentrations that appear in the image; 

Figure 16 is an amplitude spatial frequency spectrum image of the yellow-blue band of an acquired image having 
the signals of Figure 10 embedded therein; 
io Figure 17 shows geometric relationships between the peak power concentrations that appear in the amplitude 

spatial frequency spectrum image of Figure 16: 

Figure 18 shows the peak power concentrations that appear in the amplitude spatial frequency spectrum image 
of Figure 15 superimposed over the peak power concentrations that appear in the amplitude spatial frequency 
spectrum image of Figures 16 and 1 7; and 
is Figure 19 is a simplified block diagram illustrating a suitably configured machine in which the present invention 

may be used, and further illustrating the software product of the present invention and its use in conjunction with 
the machine. 

A. Adding an imperceptible signal component to an original color image. 

20 

[0019] The technique of the present invention for embedding signals in a color image has as a goal to embed the 
signals so that they are substantially imperceptible to a human viewer of the modified color image produced as the 
output of the technique. Whether the signal is perceptible is measured by a comparison operation that uses a model 
of human perception to provide a quantitative perceptual difference measurement between an original color image and 

25 a modified color image with the signal embedded therein. Perceptual difference measurements that exceed a threshold 
will be perceptible by a human viewer, and therefore cause the color image modified with the embedded signal to be 
an unacceptable degradation of the original color image. The threshold setting permits flexibility in the technique to 
accommodate the particular application that uses the technique, or to accommodate images with a particular type of 
structural or color content. Therefore, the threshold sets a limit on the acceptable perceptual difference that will be 

30 tolerated between the original color image and a color image modified with the embedded signals; the threshold may 
be changed to allow more or less perceptual difference between the original and modified color images, and therefore 
more or less degradation of the original color image. In many instances, including the illustrated embodiment described 
below, the technique can ensure that the embedded signal is imperceptible to a human viewer, and for purposes of 
this description, the embedded signal will be referred to as an imperceptible signal. It is to be understood, however, 

35 that the technique of the present invention provides a way of measuring the perceptual difference between the original 
and modified color images and to control the acceptability of this perceptual difference. 

1 . General overview and features. 

40 [0020] The technique of the present invention for embedding an imperceptible signal component in an original color 
image takes two data structures as input: the original color image (with no signal component embedded) and an image 
that contains the signal component to be added. The technique produces a version of the original color image with the 
added signal component imperceptibly embedded therein: this version of the original image is referred to herein as 
the modified image. 

45 [0021] The term "original color image" will be used to refer to the image data structure that defines the input color 
image to which the imperceptible signal is to be added Color images of the type suitable as input to the present 
invention are assumed to be two-dimensional, having an image coordinate system in which an image location is given 
by a set of image coordinates, (x, y) Each location in an image may be called a "pixel". A two-dimensional array of 
data items spatially defines the color image, each item of data providing a value, called a "pixel value", indicating the 

50 color of the image location. Each pixel value is a tuple of color space coordinates defining a color in a selected color 
space, such as RGB color space or an accepted international CIE standard color space. The original image will also 
be referred to herein as \(x,y). 

[0022] The signal component to be embedded in the original color image can itself be viewed as an image. As will 
be described in more detail below, the signal component will be embedded over several processing iterations, and so 
55 may also be referred to as S'(x,y), where / is the number of the current iteration. The characteristics of the added signal 
component vary according to the application for which the modified image carrying the signal is to be used; for some 
applications, the added signal component must be a periodic signal. In the illustrated embodiment of the technique 
that is described below, the signal added to the original color image is a sinusoidal signal component. 
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[0023] Figure 1 is a flowchart illustrating operation 200 for adding an imperceptible signal component to an original 
color image. In general terms, operation 200 handles the original color image as well as the signal component as image 
representations and adds them, in box 210, to produce modified color image, which is designated herein as V{x,y). 
[0024] To ensure that the added signal component is not perceptible in any part of the modified image by a human 

5 viewer of the image, perceptual differences between \(x,y) and V(x,y) at each image location are measured, in box 240 : 
using a comparison operation that models human perception. Comparison operation 240 produces data that indicates, 
for each image location in V(x,y), a quantitative measurement of the perceptual difference between \(x,y) and V(x,y). 
The quantitative measurement indicates whether the perceptual difference between \(x,y) and l'(x,y) is perceptible to 
a human viewer of the modified color image V(x,y). The output of comparison operation 240 is designated herein as 

10 AE{x : y) of perceptual differences. Additional details about the human perception model and a description of a model 
used in the illustrated embodiment are described below. 

[0025] The perceptual difference data is then evaluated, in box 270, in two ways. One evaluation involves whether 
any of the image locations have perceptual differences in AE(x,y) that are greater than a threshold perceptual difference 
that has been determined to be an acceptable perceptual difference measurement. Image locations with difference 

*5 measurements greater than this threshold, indicating that the added signal is perceptible to a human viewer of the 
modified color image : means that the added signal component needs to be attenuated, in a subsequent processing 
iteration, in those image locations. A second evaluation determines when sufficient iterations have been made so that 
the added signal produces no or few perceptual difference measurements over the threshold, and so is substantially 
imperceptible and processing is complete. This evaluation can be made in a number of different ways; for example, 

20 processing can be terminated when a threshold percentage of image locations in AE{x,y) have measured differences 
less than the threshold difference, or when some average (e.g., arithmetic or geometric) of the AE values is below a 
threshold. Both of these thresholds can be adjusted to suit the specific application in which the modified color image 
is to be used, or to accommodate the particular spectral or structural attributes of the original color image. 
[0026] As box 270 in Figure 1 indicates : when both evaluations show that the modified color image has a sufficient 

2S number of image locations with difference measurements below the threshold acceptable difference, operation 200 
for producing the modified color image is complete, and processing stops. The final version of modified color image 
290 has an added signal component that produces an acceptable perceptual difference between the original color 
image and the modified color image; in most cases, the added signal component is imperceptible to human viewers 
of the image 

30 [0027] When image locations in l'(x,y) have perceptual differences in AE(x,y) that are greater than the threshold 
acceptable difference measurement, a new signal image is produced, represented in box 280 of Figure 1 , in which the 
- signal of the previous signal image is attenuated in image areas where there were unacceptable difference measure- 
ments. This new signal image is produced using the signal image S(x,y) of the previous iteration, to form the signal 
image, designated as S> +1 (x,y) herein, for the next (/+ 1) iteration. Processing then returns to box 210 where modified 

35 signal S* +1 (x,y) is added to the original image. 

[0028] Figures 2 - 6 illustrate operation 200 in terms of representative schematic image representations. Original 
color image 300 of Figure 2 shows image regions 302, 304, 306, 308 and 310, each shown in a different color, as 
represented by different cross-hatching patterns. Figure 3 shows image 320 representative of an image of the signal 
component to be added to original image 300. Note that image 320 is shown with a cross-hatching pattern that is 

^0 intended to be representative of any type of signal component, in contrast to image 480 of Figure 10 which is an image 
of a specific signal component that is the subject of an illustrated embodiment described in more detail below. 
[0029] Figure 4 shows image 330 produced by operation 210 of Figure 1 which is representative of the image pro- 
duced by adding image 320 to image 300. In effect, the image produced as a result of the process of embedding signals 
in the original color image according to the present invention is defined to be two image data structures combined and 

45 rendered in registered alignment. In image 330, image regions 332, 334 and 340 are shown having the same colors 
(i.e., with the same cross-hatching) as in image 300 of Figure 2, while image regions 336 and 338 are shown with the 
cross-hatching color of signal image 320 superimposed over the cross-hatching colors of image regions 306 and 308. 
Image 330 schematically illustrates that signal component 320 is perceptible in image regions 336 and 338 but not in 
image regions 332, 334 and 340. Figure 5 shows image 350 representative of perceptual differences, AE(x,y), produced 

50 by operation 240 which measures the perceptual differences between image 300 and image 330 using the human 
perception model. The quantitative perceptual difference AE in image regions 358 and 356 are above the acceptable 
perceptual difference threshold. 

[0030] Figure 6 shows image 360 representative of the signal image S* +1 (x,y) produced by operation 280 of Figure 
2. Signal image 360 is to be embedded in the original image in the next iteration of processing, and shows, by way of 
55 the different cross-hatching pattern in region 364, that the signal component is attenuated in the image regions corre- 
sponding to the image regions in image 330 where the added signal component produces an unacceptably large 
perceptual difference and is therefore perceptible to a human viewer of image 330. 
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2. Implementation details of an illustrated embodiment. 

[0031] An illustrated embodiment of the operation 200 of the present invention adds to an original color image one 
or more sinusoidal grids as the signal component where the sinusoids have predetermined spatial directions and 

5 spatial frequencies. The illustrated embodiment is shown in flowchart 400 in Figure 7. Original color image 404 is 
converted, in box 406, to an opponent color representation which produces three color separations: luminance (black 
and white) : red-green (RG) and yellow-blue (YB). This conversion is a linear transformation, typically from RGB color 
space. A sum of sinusoidal signals 402, S(x,y), is embedded into the yellow-blue (YB) color band in box 408. The 
illustrated embodiment chooses high spatial frequencies where human sensitivity is low and where most images contain 

10 relatively tittle power. The three color separations including the YB band with the embedded signal form an interim 
modified image 414 which is then converted to an input representation suitable for input to the S-CIELAB perceptual 
difference measurement operation 450. In particular, interim color image 414 is converted, in box 409, to a represen- 
tation that approximates the absorptions of the long, middle and short (LMS) wavelength sensitive cone photoreceptors 
in the retina of the human eye. This LMS representation, along with the LMS representation of original image 404, also 

is converted in box 409, are input to comparison operation 450 to detect and measure perceptual differences between 
the two images. When the perceptual difference measurements exceed the thresholds, the sinusoidal grid 402 must 
be attenuated in those areas that produce perceptible differences, and again added to the YB band of original image 
404. Thus, the signal is embedded in an iterative manner. The product of sinusoidal grid image 402 and an attenuation 
image 41 8 produced in operation 420 is added to original image 404. If the perceptual difference measurements again 

20 exceed the thresholds, another iteration is needed; a new attenuation image 418 is produced in box 420, which is again 
multiplied with sinusoidal grid 402 produced in the last iteration, to form the new sinusoidal grid to be added to the YB 
band of the original image 404 in the next iteration. The iterative processing is represented by line 410. 
[0032] Note that, while the technique of the present invention refers to embedding signals in color images : the signal 
component is not restricted to being embedded in only the color (i.e., RG or YB) bands of the opponent-color image 

25 representation. While the technique of the illustrated embodiment embeds the signal component in the YB band of the 
original image, the signal may be embedded in any one or more of the BW, RG or YB bands of the original image The 
perceptibility of the embedded signal in a particular one of the opponent color representations may depend on the 
spatial or color structure of the particular image. Note also that other types of image separations are also possible, 
such as. for example, HSV (hue, saturation and value) and L*a*b*. 

30 [0033] The illustrated embodiment uses a human perception model known as S-CIELAB to make the perceptual 
difference comparisons between the original image and each iterative, interim version 414 of the modified image, 
producing the perceptual difference measurement image 416, AE(x,y). The functional and operational characteristics 
of S-CIELAB are described below. Note that the conversion of an RGB image to its LMS representation required in 
box 409 of Figure 7 is a standard operation in the literature of human perception, and is provided as a utility with the 

35 S-CIELAB source code. 

[0034] The technique of the present invention, including the process of embedding the sinusoidal grid and measuring 
the perceptual differences, as well as a decoding process to determine the location of the sinusoidal grid in an acquired 
modified color image, have been implemented in the programming language known as Matlab; these implementations 
execute on a suitably configured general purpose computer, such as machine 1 00 schematically represented in Figure 

40 t9 ) and described below. 

i. Human color perception model. 

[0035] As noted earlier, the purpose of comparison operation 240 of Figure 1 is to measure the perceptual differences 
45 between the original color image and the current version of the modified color image to determine if the perceptual 
differences exceed the threshold of acceptable perceptual differences and are consequently perceptible by a human 
viewer of the image. Many models of human color perception have been proposed in the color research literature; 
models that are suitable for use with the present invention must provide the functional ability to measure a color ap- 
pearance difference between two colors that represents the perceptual difference that a human observer would per- 
50 ceive between the two colors, and to express the measured difference in a quantitative range, or scale, of small to 
large perceived color appearance differences that accurately corresponds to actual human color vision. The technique 
of the present invention must be able to select a threshold perceptual difference from the scale of quantitative differ- 
ences produced as a result of using the model that represents a perceptual difference that is "acceptable", either 
because the threshold represents a perceptual difference that falls below the threshold of human perception, or because 
55 perceptual differences at or above the threshold are perceptible to a degree that is acceptable for the particular appli- 
cation. In addition, a model that is suitable for use with the present invention wilt premise its color measurements on 
the actual physical responses of the human eye to color, and will generally allow for specification of the white and black 
points in the image and/or an illuminant used to produce the image. CIE (Commission Internationale de I'Eclairage) 
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standards for defining and measuring colors form the basis for many such models of human color perception. Note 
that how well a particular human color perception model accurately represents and measures perceptual differences 
is likely to influence the operation of the present invention; the more accurate information about perceptual differences 
that comparison operation 240 of Figure 1 provides, the more likely the technique can ensure that the embedded signal 
5 is imperceptible. 

[0036] In the illustrated embodiment, the present invention makes use of a human color perception model known as 
S-CIELAB. The S-CIELAB model was developed by X. Zhang and B. Wandell, and is based on the color vision research 
of A. Poirson and B. Wandell. Additional information about S-CIELAB and Poirson and Wandell's research beyond that 
described below is available in: A. Poirson and B. Wandell, "The appearance of colored patterns: pattern-color sepa- 

io rability", Journal of the Optical Society of America A, 10:2458-2471, 1993; A. Poirson and B. Wandell, "Pattern-color 
separable pathways predict sensitivity to single colored patterns", Vision Research, 36:515-526, 1996: X. Zhang and 
B. Wandell. "A spatial extension of CIELAB for digital color image reproduction", Proceedings of the Society of Infor- 
mation Display 96 Digest, pp. 731 -734 : San Diego, 1996. Each of these publications provides references to color 
science and human vision texts that may be consulted as necessary. Source code for an S-CIELAB implementation 

is is available for downloading from Brian Wandell's web site (URL at filing date: http://white.stanford.edu/html/wandell. 
html). 

[0037] S-CIELAB is an extension of the well-known CIELAB standard that attempts to define a perceptually uniform 
color space in which color reproduction errors can be accurately measured. S-CIELAB was designed to measure the 
perceptual difference errors between an original color image and its color reproduction on a particular device. S-CIELAB 

20 attempts to account for the dependence of color appearance on the spatial structure of an image, both globally and 
locally. CIELAB was based on psychophysical studies of color discriminability with large uniform color patches. More 
recent studies have shown that color appearance also depends on the local structure of the image. Measurements of 
color appearance with fine-scale patterns (e.g. , color square-wave gratings of different spatial frequencies) differ from 
those with large uniform color patches. Poirson and Wandell propose a pattern-color separable model for early human 

2S color processing to explain these human color vision responses. Zhang and Wandell state that applying CIELAB to 
predict local color reproduction errors in patterned images does not give satisfactory results, and cite as an example 
the comparison of a continuous-tone color image with a halftone version of the image. A point-by-point comparison of 
the CIELAB error between these images produces large errors at most image points. However, because the halftone 
patterns vary rapidly, these differences are blurred by the eye, and the reproduction may still preserve the appearance 

30 of the original. 

[0038] Figure 8 illustrates the general processing of S-CIELAB comparison operation 450 as implemented in the 
illustrated embodiment using the S-CIELAB color perception model. In the S-CIELAB pattern-color separable model, 
the color signals of both the original and the current version of the modified color images are input in an LMS repre- 
sentation. The color signals in this representation are then transformed, in box 454, from LMS coordinates into an 

35 opponent-color representation having three bands (i.e., spatial images), referred to as luminance (BW), red-green 
(RG), and yellow-blue (YB). Each opponent-color image is then convolved, in box 456, with a kernel whose shape is 
determined by the visual spatial sensitivity to that color dimension: the area under each of these kernels integrates to 
one. This amounts to smoothing each opponent-colors image with a lowpass filter, where the RG and YB bands have 
lower frequency cutoffs. The calculation is pattern-color separable because the color transformation does not depend 

40 on the image's spatial pattern, and the spatial convolution does not depend on the image's color. 

[0039] Finally, the filtered representation of the color signals at each image location are transformed into CIE XYZ- 
tristimulus coordinates, in box 458. The CIE XYZtristimulus representation is then transformed, in box 460, into CIELAB 
color space using the CIELAB formulae. The conventional CIELAB color difference lormula is then used to measure 
the perceptual difference at each image pixel between a color in the original color image and a color in the current 

45 version of the modified color image. These perceptual differences are denoted by AE, where values greater than 1 
imply that the perceptual differences are detectable by humans. A perceptual difference value of less than 1 implies 
that a difference is not discriminate. The perceptual difference image, AE(x,y) : is then returned to operation 420. The 
S-CIELAB difference measure reflects both spatial and color sensitivity, and, according to Zhang and Wandell, it equals 
the conventional CIELAB over uniform regions of the image. 

so [0040] Use of the S-CIELAB color perception model requires some calibration parameters. S-CIELAB requires the 
mapping between the digital representation of the color image and the spectral reflectance of the image when printed 
or displayed. Thus, a detailed model of the printer or display device and the spectral distribution of the illuminant need 
to be specified to the model. In the illustrated embodiment of the present invention, the S-CIELAB model assumes that 
the image is displayed on a CRT monitor with a conventional set of RGB phosphors The S-CIELAB model also requires 

55 the specification of a white point. In the illustrated embodiment of the present invention, the S-CIELAB model assumes, 
for the digital images being processed, that the white point is the color tuple at that pixel with the highest luminance in 
the original color image; that is, the largest Y component, in CIE XYZ-tristimulus coordinates, of all of the image pixels 
is designated as the white point. For images printed on white paper, the white point is normally specified as the reflect- 
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ance of the paper alone. 

[0041] Experience with the S-CIELAB model shows that, in practice, perceptual differences with AE< 3 in the modified 
color image with the embedded sinusoid signal component are not readily apparent. However, when, as in the case 
of the illustrated embodiment, there is incomplete knowledge of the proper calibration parameters, it is important to 
5 keep the AE values small (e.g., < 1 ) to be conservative. As noted earlier, the threshold can be changed to accept more 
image degradation in the final modified color image. 

ii. Producing the attenuated signal component. 

to [0042] As noted earlier, in the illustrated embodiment of the present invention shown in Figure 7, perceptual difference 
image 41 6, AE{x,y) represents the perceptual differences between original image 404, l(x.y), and the current, or interim, 
modified image 41 A, l'(x,y). The signal to be embedded is attenuated wherever the perceptual differences between 
the original image and the modified image is significant, i.e., where AE(x,y) > 1. The flowchart of Figure 9 illustrates 
operation 420 for producing the signal component image, S^{x,y), for the next iteration, according to the illustrated 

is embodiment. At each iteration ii there are a sufficient number of image locations where AE(x,y) > t, where / is some 
threshold, then AE{x,y) is converted into an image of attenuation factors a(x,y), in box 422, also called attenuation 
image 418 (Figure 7), as follows: 



20 



35 



40 



a{x,y) = (2/2 + max[0, AE{x } y) - t]f (1) 



[0043] The threshold t thus represents the threshold below which no A£(x,y) will cause attenuation of the signal 
image. In the illustrated embodiment, the threshold used is 1. 

[0044] The image of attenuation factors, a(x,y), is then blurred with a lowpass filter, in box 426. The blurring of a(x, 
25 y) is important because S-CIELAB measures perceptual differences in local neighborhoods, due to the spatial blurring 
in each opponent-color band. The image of attenuation factors, a(x,y) : is blurred with the same lowpass filter that S- 
CIELAB uses for the YB band. The blurring also keeps the bandwidth of the embedded signal relatively small. This 
keeps its frequency content away from low frequencies to which humans are very sensitive. Smaller bandwidths also 
keep the power of the sinusoids more concentrated and easier to detect. 
30 [0045] The blurred attenuation image 418 is then multiplied, in box 428, with the signal that was embedded in the 
previous iteration. Operation 400 of the illustrated embodiment (Figure 7) repeats until 99.5% of the pixels haveA£:(x, 
y) □ t, which has been found to take 3 to 5 iterations. The final embedded signal is given by 



G{x, y ) = (IT aiUy)) s&y) (2) 
1=1 



where nis the number of iterations, S(x,y) is the sum of sinusoids, and a(x,y) is the blurred attenuation image at iteration 
/. The amplitude of the embedded signal can only decrease at each iteration, which ensures convergence so that AE 
(x,y) □ 1 almost everywhere. Equation (1), which converts AE values into attenuation factors, was chosen empirically 
to obtain a reasonable rate of convergence. Attenuation of the sinusoidal signals must not fall below a limit below which 

45 the signals would not be detectable by an image acquisition device such as a scanner. 

[0046] Selection of the number of sinusoids to embed and the sinusoid spatial directions and spatial wavelengths 
may be influenced by the type of application in which the modified image, or a version of the modified image, will be 
used. In the illustrated embodiment of the present invention, four sinusoids were embedded into a color image, forming 
a sinusoidal grid over the image which, when decoded 1rom an acquired image, may be used to determine certain 

50 geometric relationships between the acquired image and the modified image The original image had R, G and B values 
at each pixel each expressed in the range of 0-255 (i.e. , 8 bits per color band per pixel). After conversion to the opponent 
color representation, the YB (yellow-blue) channel typically has values between -8 and 8. Four sinusoids are added 
with amplitudes of approximately .5 pixels each. Therefore, the typical amplitude of the embedded signal in this example 
was equal to 2. 

55 [0047] A schematic representation of these four sinusoidal signals that form the signal image S(x,y) is illustrated in 
spatial image 480 of Figure 10. Image 480 shows pairs of equally-spaced diagonally positioned lines extending across 
image 480 in orthogonal directions, forming an image of a grid. The lines mark the zero crossing positions of sinusoids 
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having different spatial frequencies and positioned in different spatial directions in image 480. It can be seen trom the 
relative positrons of the pairs of diagonal lines that the four sinusoids have orientations of ±45° from a coordinate axis 
(not shown) in the x direction of image 480. Figure 13 shows diagonal cross-section 482 of image 480, showing the 
relationship of a pair of sinusoids; dashed line 481 marks the horizontal sinusoidal axis. In the illustrated embodiment, 
the wavelengths of the sinusoids are 3.8 and 2.7 pixels, represented by wavelengths 486 and 488. S-CIELAB was 
calibrated so that the embedded signal would be invisible when viewed from 10 inches, at which the frequencies are 
7 and 10 cycles/degree. In actual experience, the process was sufficiently robust that the embedded signal remained 
imperceptible even when viewed more closely than 10 inches. 

[0048] In another application, the desired result of decoding the sinusoidal image from an acquired image version 
of a modified image may be to facilitate establishing absolute spatial coordinates in the acquired image. This may be 
accomplished by adding two or more sinusoids to an original image that have specific characteristics. In particular, the 
added sinusoids should each have a controlled, predetermined phase, and the phase of the first sinusoid should be 
different from that of the second sinusoid. The total combined signal formed by the two sinusoids should have a period 
greater than that of the image extent; so, for example, if the sinusoids have a horizontal spatial direction, the total 
combined period should be greater than the horizontal dimension (e.g., width) of the image. The relative phases of the 
sinusoids combined with the fact that the total period of the signal is greater than the image extent allow for designating 
the location where the two (or more) sinusoids simultaneously cross their horizontal axis to be a unique spatial image 
origin in the acquired image. Figure 15 shows sinusoidal image 490 which shows sinusoids 492 and 494. Each of these 
sinusoids has a phase different from the other and a period greater than the width of the image; they simultaneously 
cross their common axis at zero only at point 496. 

B. Applications for color images with embedded signals. 

1. Overview of application processing. 

[0049] In this discussion, the term "acquired image" is used to mean an image that has been digitally captured by a 
device from a visual rendering of a modified image produced as the output of either operation 200 (Figure t ) or operation 
400 (Figure 7). Since the modified image is a version of the original image with the periodic signals embedded therein, 
the modified image acts as a surrogate for the original image, and applications which make use of the signals embedded 

30 in the acquired image do so in order to perform operations or computations with respect to the original image. The 
modified image, and acquired images that are versions of the modified image, in effect carry with them information 
about the original image without unacceptably altering the aesthetic appearance of the original image. 
[0050] Typically, an acquired image is produced by a scanning device that optically senses the colors (black and 
white or chromatic) of a printed rendering of the modified image and converts the captured signals into a digital rep- 

35 resentation of the image. However, other types of digital capture devices, such as a digital camera, might also be the 
signal source for the acquired image. 

[0051] Figure 13 provides a general flowchart 500 of a class of applications in which an image modified according 
to the technique illustrated in Figures 1 or 7 may be used. In order to make clear the distinction between acquired 
image 504 and modified image 501 in this discussion, the image acquisition process is schematically represented in 
^0 Figure 13 as operation 502: the dotted lines surrounding image acquisition operation 502 indicate that it is a process 
distinct and separate from the class of applications illustrated in flowchart 500. 

[0052] The class of applications illustrated in flowchart 500 uses signals imperceptibly embedded in an acquired 
image to determine geometric differences between the acquired image and the original image, without using the actual 
original image during the decoding process. The signal embedded in the modified image is a periodic signal; information 

45 decoded from the acquired image about the spatial frequencies of the embedded periodic signal components along 
with predetermined geometric information about the embedded signals are used to determine a geometric relationship 
between the embedded signals in the acquired image and the original embedded signal image. This geometric rela- 
tionship between the signals may, in turn, be used to compute and automatically match the geometric properties of 
the acquired image to those of the original image without requiring the input of either the modified or the original image. 

so Such geometric properties include orientation, scale, and linear deformations. Additionally, the phases of the embedded 
signals in the acquired image can be used to establish absolute local coordinates in the acquired image from which to 
determine image positions. 

[0053] Flowchart 500 illustrating this class of applications proceeds as follows. An acquired image 504 and geometric 
constraint information 508 about the embedded signals are input to decoding operation 520, which determines the 
55 spatial frequencies of the embedded signals in acquired image 504. The output of the decoding operation includes a 
linear mapping 558 of locations of local peak power concentrations between acquired image 504 and modified image 
290. This linear mapping information is input to operation 560 which computes certain geometric relationships between 
acquired image 504 and modified color image 290. Once these geometric relationships are known, acquired image 
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504 may be geometrically transformed, in box 566 to match the geometric properties of modified color image 290, 
producing geometrically modified acquired image 570 

2. Decoding an acquired image to locate the added signal component. 

5 

[0054] Figure 14 provides a general flowchart of an illustrated embodiment of decoding operation 520. The goal of 
decoding operation 520 is to locate and identily the spatial frequencies of the sinusoidal signals embedded in an 
acquired image in the manner described above with respect to the illustrated embodiment of the signal embedding 
technique shown in Figures 7 and 10. Image 480 of Figure 10 including the set of four (4) sinusoids has a distinctive 

10 pattern of local peak power concentrations indicating the locations of the spatial frequencies of the sinusoids; this 
pattern can be seen in Figure 15. Figure 15 shows image 580 of an enhanced log amplitude spatial frequency spectrum 
of the signal component image of Figure 10. Image 580 has an image power spectrum origin at point 581 and has 
been thresholded to show the local peak power concentrations (hereafter simply referred to as peaks) represented at 
points 582, 583, 584, 585, 586, 587, 588 and 589. It can be seen that these peaks occur in two sets of four each and 

is have geometric relationships to each other with respect to coordinate axis 578 having an origin at position 581. Spe- 
cifically, the peaks occur at distances 576 and 577 from position 531 (as measured along dashed line 574) and at an 
angle 575 above and below coordinate axis 578. It was previously noted in the discussion of Figure 1 0 that the sinusoids 
have orientations of ±45° and wavelengths of 3.8 and 2.7 pixels. Image 580 shows the peaks corresponding to these 
specific characteristics of this signal. In more general terms, an image composed of a set of n sinusoids has the geo- 

20 metric properties of angular orientation differences and radial frequency ratios from a point assigned as an origin of 
the power spectrum image. These geometric properties of the embedded signal are referred to herein as "geometric 
constraints", and are the only apriori information that is needed to decode the spatial frequencies from an acquired 
image in decoding operation 520. The source of the geometric constraint data is unimportant; it may be made available 
as separate input data from another signal source; or it may be included in the acquired image in some manner. 

25 [0055] Returning to Figure 14, acquired image data structure 504 is first converted into a plurality of opponent-color 
representations; this is accomplished in the same manner as the original and interim color images were converted into 
opponent-color representations, as described earlier in the discussion of illustrated embodiment 400 of Figures 7 and 
8, in boxes 409 and 454. An FFT operation 524 is applied to the YB color separation to compute the amplitude spatial 
frequency spectrum of the YB band; this spectrum is then divided by a blurred version of itself, in box 525, to enhance 

30 the local peak power concentrations. Then, all power peak concentrations in the YB band of the acquired image are 
located, in box 526, and included as candidate spatial frequencies to be tested. 

[0056] Figure 16 shows image 590 of an enhanced log amplitude spatial frequency spectrum of the YB band of 
acquired image 504 having the sinusoid image 480 of Figure 10 embedded therein. Image 590 has an image origin at 
point 591 and been thresholded to show all local peaks. With reference again to Figure 14, sets of 4 local peaks are 

35 compared against the geometric constraints 508 that are shown in Figure 1 5, specifically the distance ratio of distances 
576 and 577 and angle 575. Those sets of four frequencies that violate constraints 508 are discarded, in box 528. Of 
the remaining candidate frequencies, the best linear mapping is found, in box 530, between the locations of the local 
peak power concentrations in the acquired image and the locations of the local peak power concentrations in the 
sinusoidal image, using the geometric constraints 508. A least squares technique may be used to produce the best 

40 mapping. Box 530 produces linear mapping 558 as its output, and decoding operation 520 is complete. 

3. Automatic image alignment and sealing. 

[0057] As noted in the discussion accompanying Figure 1 3, in this class of applications tor an image having embed- 
45 ded signals, the linear mapping information is input to an operation that computes geometric relationships between an 
acquired image and an original image. With respect to the illustrated embodiment of the amplitude spatial frequency 
spectrum of the acquired image shown in Figure 16, processing in box 530 (Figure 14) shows that the best match of 
a set of 4 spatial frequencies from the remaining candidate frequencies are found at points 592, 593, 594, 595, 596, 
597, 598 and 599. As shown in Figure 17, once these peaks are identified, the angle 614 of orientation from an image 
so coordinate axis 610 and the distances 618 and 620 of the peaks away from point 591 can be computed. It can be seen 
from a visual comparison between Figures 15 and 17 of the position and orientation of these peaks that the peaks in 
the acquired image as represented in Figure 17 have geometric properties that differ from those of the peaks in the 
sinusoidal image as represented in Figure 1 5. Figure 1 8 illustrates these geometric differences. In Figure 18, amplitude 
spatial frequency spectrum image 580 of the sinusoidal image is superimposed on amplitude spatial frequency spee- 
ds trum 590 of the acquired image and aligned at their designated origin points. The local peaks in image 580 are shown 
with dashed lines connecting them, while the local peaks in image 590 are shown with solid lines connecting them. It 
can be seen that the peaks are not aligned; the peak power concentrations in the amplitude spatial frequency spectrum 
590 of the acquired image are offset by angle 650, and there are differences in distances of the peaks from the origin 
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by distances 652 and 654. Angle 650 difference represents a rotation of acquired image 504, which may have resulted 
from misplacement of the visual rendering of the modified image during the image acquisition process. Distance dif- 
ferences 652 and 654 represent a scaling (enlargement) of acquired image 504, which may have occurred, for example, 
during the rendering of the modified color image. Linear mapping 558, produced as the output of decoding operation 
5 520, is used to compute these geometric differences. Once these geometric properties are computed, acquired image 
504 may be automatically reduced in scale and rotated to match the geometric properties of the original image. The 
class of applications illustrated in the flowchart of Figure 1 3 include all such linear geometric transformations between 
an acquired image and an original image having embedded periodic signals. 

10 4. Other applications for color images having embedded signals. 

[0058] As discussed in the background section, much research and development is concerned with imperceptibly 
embedding information in images. The embedding technique of the illustrated embodiment may be used in conjunction 
with any other suitable information embedding technique. As noted earlier the sinusoids of Figure 10 when embedded 

15 in an acquired image produce a grid across the image. The decoding operation locates and identifies the spatial fre- 
quencies; this information along with the geometric constraint information of the original image enables the acquired 
image to be geometrically matched to the original image. Decoding ot the phase of the embedded sinusoids enables 
an absolute coordinate frame to be assigned to the acquired images, provided that the sinusoids are embedded as 
described in the discussion accompanying Figure 11 . Knowledge of the location of the sinusoidal grid in the acquired 

20 image permits the grid to function as a landmark for other information that may be imperceptibly embedded in the 
acquired image. The ability to orient and locate other embedded information is another important application of the 
embedding technique described herein. 

[0059] Another feature of the present invention is that the technique for decoding the sinusoidal grid can be straight- 
forwardly extended to determine whether there are locally smooth nonlinear deformations of the modified image, such 
2S as those that occur when the image has not been placed flat on the platen of a scanner. Detection of these deformations 
produces local geometric inhomogeneities of the sinusoidal grid; the ability to map the modified image to the original 
image allows automatic correction of the deformed portion(s) of the image. 

D. The Machine and Software Product of the Invention. 

30 

[0060] Figure 19 is a block diagram of a generalized, processor-controlled machine 100; the present invention may 
be used in any machine having the common components, characteristics, and configuration of machine 100, and is 
not inherently related to any particular processor, machine, system or other apparatus. The machine or system may 
be specially constructed and optimized for the purpose of carrying out the invention, or it may comprise a general 

35 purpose computer selectively activated or reconfigured by a computer program stored in the computer, or it may be a 
combination of a general purpose computer and auxiliary special purpose hardware. When a machine such as machine 
100 is suitably programmed to embody the present invention, the machine is not a standard or known configuration. 
[0061] Machine 100 may be configured to perform any one of operation 200 of Figure 1 to embed data in a color 
image, operation 520 of Figure 14 to decode an embedded signal from an acquired image, or operation 500 of Figure 

40 1 3 to compute the geometric relationship between an acquired and a modified color image, or it may be configured to 
perform all of these operations. An input color image, such as original color image color image 300 of Figure 2, color 
image 404 of Figure 7, or acquired image 504 of Figure 1 3, is provided from signal source 1 58. Signal source 1 58 may 
be a memory device, a communications channel, a data bus : another processor performing an operation, or any other 
suitable source of color image signals. For example, signal source 158 may be an image capture device, such as a 

45 scanning device, a digital camera, or an interface device that produces a digital image definition data structure from 
another type of image signal. An input color image provided by signal source 156 is forwarded via input circuitry 156 
to processor 140 and may be stored in data memory 114. Machine 100 may : but need not, include a conventional 
display device (not shown) capable of presenting images, such as a cathode ray tube, a liquid crystal display (LCD) 
device, a printing device, or any other device suitable for presenting images. 

so [0062] Processor 140 operates by accessing program memory 110 to retrieve instructions, which it then executes. 
When machine 100 is configured to perform signal embedding operation 200, program memory 110 includes signal 
embedding instructions 200 that implement the functions shown in flowchart 200 of Figure 1 Program memory 110 
includes instructions for the subroutines needed to produce the modified color image having the embedded signals 
according to signal embedding instructions 200; these include perceptual difference measurement instructions 240 

55 and attenuation image production instructions 280. Similarly, when machine 100 is configured to perform embedded 
signal decoding operation 520, program memory 110 includes decoding instructions 520 that implement the functions 
shown in the flowchart of Figure 14, and also includes instructions for the subroutines needed to decode the spatial 
frequencies from an acquired image. When machine 100 is configured to perform automatic image scaling and align- 
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ment, or another type of application, on an acquired image, program memory 110 includes application instructions 500 
that implement the functions shown in the flowchart of Figure 13, including instructions for the subroutines shown in 
Figure 13 for transforming the geometric properties of the acquired image to match those of the original color image. 
[0063] During execution of the instructions, processor 140 may access data memory 114 to obtain or store data 

5 necessary for performing its operations. For example, when machine 100 is configured to perform operation 200 for 
embedding a signal in an original color image, data memory 114 stores the image data structure 404 defining the 
original color image as well as the image data structures defining the opponent-colors representations 412 of the 
original and interim versions of the modified color images, and the final modified color image 424 with the embedded 
signals. Data memory 114 also stores signal component data 408 to be added to original color image 404, and data 

10 416 defining perceptual difference measurements AE(x,y). Data memory 114 also stores various other miscellaneous 
data. Similarly, when machine 100 is configured to perform decoding operation 520, data memory 114 stores geometric 
constraints 508, candidate spatial frequencies, and linear mapping 558 of locations of local peak power concentrations 
between acquired and modified image spatial frequency power spectra. 

[0064] The actual manner in which the physical components of machine 100 are connected may vary, and may 
is include hardwired physical connections between some or all of the components, as well as connections over wired or 
wireless communications facilities, such as through remote or local communications networks and infrared and radio 
connections. Program memory 110 or data memory 114, for example, may include memory that is physically connected 
to processor 140 as local memory, or that is remotely accessible to processor 140 by means o1 a wired or wireless 
communications facility (not shown). 
20 [0065] Figure 19 also shows software product 160, an article of manufacture that can be used in a machine that 
includes components like those shown in machine 100. Software product 160 includes data storage medium 170 that 
can be accessed by storage medium access circuitry 150. Data storage medium 170 stores instructions for executing 
one or more of operation 200 of Figure 1, operation 500 of Figure 13, or operation 520 of Figure 14, and may also 
include instructions for performing these operations according to the illustrated embodiments illustrated in Figures 7, 
2S 8 and 9. 

[0066] Software product 160 may be commercially available to a consumer in the form of a shrink-wrap package 
that includes data storage medium 170 and appropriate documentation describing the product. In that case, a data 
storage medium is a physical medium that stores instruction data. Examples of data storage media include magnetic 
media such as floppy disks, diskettes and PC cards (also known as PCMCIA memory cards), optical media such as 

30 CD-ROMs, and semiconductor media such as semiconductor ROMs and RAMs As used herein, "storage medium" 
covers one or more distinct units of a medium that together store a body ol data. For example, a set of disks storing 
a single body of data would be a storage medium. "Storage medium access circuitry" is circuitry that can access data 
on a data storage medium. Storage medium access circuitry 150 may be contained in a distinct physical device into 
which data storage medium 170 is inserted in order for the storage medium access circuitry to access the data stored 

35 thereon. Examples of storage medium access devices include disk drives and CD-ROM readers. These may be phys- 
ically separate devices from machine 100, or enclosed as part of a housing of machine 100 that includes other com- 
ponents. 

[0067] Storage medium access circuitry 150 may also be incorporated as part of the functionality of machine 100, 
such as when storage medium access circuitry includes communications access software and circuitry in order to 

40 access the instruction data on data storage medium 1 70 when data storage medium 1 70 is stored as part of a remotely- 
located storage device, such as a server. Software product 160 may be commercially or otherwise available to a user 
in the form of a data stream indicating instruction data for performing the method of the present invention that is trans- 
mitted to the user over a communications facility from the remotely-located storage device. In the latter case, article 
160 is embodied in physical form as signals stored on the remotely-located storage device; the user purchases or 

45 accesses a copy of the contents of data storage medium 170 containing instructions for performing the present inven- 
tion, but typically does not purchase or acquire any rights in the actual remotely-located storage device. When software 
product 160 is provided in the form of a data stream transmitted to the user over a communications facility from the 
remotely-located storage device, instruction data stored on data storage medium 170 is accessible using storage me- 
dium access circuitry 150. Alternatively, a data stream transmitted to the user over a communications facility from the 

so remotely-located storage device may be stored in some suitable local memory device of machine 100, which might be 
program memory 1 1 0, or a data storage medium locally accessible to machine 1 00 (not shown), which would then also 
be accessible using storage medium access circuitry 150. 

[0068] Figure 19 shows data storage medium 170 configured for storing instruction data for performing operations 
200 (Figure 1), 500 (Figure 13) and 520 (Figure 14). This instruction data is provided to processor 140 for execution 
55 when the technique for embedding imperceptible signal data in a color image is to be used. The stored data includes 
binary data signal embedding instructions 162 and 164; when these instructions are provided to processor 140, and 
processor 140 executes them, the machine is operated to perform the operations for adding imperceptible signal data 
to an original color image, as represented in the flowchart of Figure 1 . For example, as shown in Figure 19 : the stored 
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data include perceptual difference measurement instructions 164; when these instructions are provided to processor 
140, and processor 1 40 executes them, the machine is operated to perform the operations for measuring perceptual 
differences between the original color image and the current version of the modified color image using a human per- 
ception model, as represented in box 240 of Figure t, or in the flowchart of Figure 8. The stored data shown also 
s include attenuation image production instructions 162; when these instructions are provided to processor 140 : and 
processor 140 executes them, the machine is operated to perform the operations for producing an attenuated signal 
component for the next iteration of the signal embedding operation, as represented in box 280 of Figure 1, or in the 
flowchart of Figure 9. 

[0069] Also shown in Figure 19, the data stored on data storage medium 170 further includes instruction data 166 
to indicating decoding instructions for decoding the signal component embedded in a modified color image; when these 
instructions are provided to processor 140, and processor 140 executes them, the machine is operated to perform 
decoding operation 520, as represented in the flowchart of Figure 14. The data stored on data storage medium 170 
further includes instruction data 168 indicating geometric transformation instructions for automatically transforming the 
geometric attributes, such as skew and scale, of an acquired image to match those of its corresponding original color 
is image; when these instructions are provided to processor 140, and processor 140 executes them, the machine is 
operated to perform automatic spatial distortion correction, deskewing and scaling of an acquired image, as one of the 
class of applications represented in operation 500 of Figure 13 



20 Claims 

1. A method for operating a machine to produce a version of an original color image data structure having perceptually 
acceptable signal differences embedded therein; the machine including a processor and a memory device for 
storing data; the data stored in the memory device including instruction data the processor executes to operate 

25 the machine; the processor being connected to the memory device for accessing and executing the instruction 

data stored therein; the method comprising: 

operating the processor to obtain an original color image data structure defining an original color image in- 
cluding a plurality of original image locations each indicating an original color; 

30 operating the processor to obtain a signal; 

operating the processor to add the signal to the original color image data structure to produce an interim color 
image data structure; the interim color image data structure being a version of at least a portion of the original 
color image data structure and including a plurality of interim image locations each indicating an interim color; 
the signal being added to at least one of the plurality of interim image locations; an original image location in 

35 the original color image data structure and an interim image location in the interim color data structure indicating 

the same image location in each image and being referred to collectively as respectively paired image loca- 
tions; each one of the respectively paired image locations indicating an image color: 

operating the processor to perform a comparison operation to produce a quantitative perceptual difference 
measurement at each of the respectively paired image locations; the comparison operation using a model of 
40 human perception to quantitatively represent the perceptual differences as perceived by a human viewer of 

the original and interim color images; 

operating the processor to attenuate the signal to produce a non-zero modified signal when the difference 
measurements between the respectively paired image locations exceed a threshold indicating an acceptable 
perceptual difference between the original and interim color images; and 
45 operating the processor to perform an iteration operation iterating the adding, comparison and attenuating 

operations using the modified signal until the perceptual difference measurements between the respectively 
paired image locations do not exceed the threshold; the iteration operation producing a modified color image 
data structure indicating a perceptually acceptable signal difference between the original and modified color 
images. 

so 

2. A method according to claim 1, further including, prior to adding the signal to the original color image structure, 
operating the processor to separate the original color image data structure into a plurality of individual color rep- 
resentations; and wherein the signal is added to a selected one of the individual color representations. 

55 3. A method according to claim 1 or claim 2, wherein the model of human perception models human perception of 
an image color with respect to a local spatial structure in which the image color is located in an image; the quantative 
perceptual difference measurement indicating a perceptual difference that depends on at least one of the local 
spatial structure or color of the image. 
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4. A method according to any of the preceding claims, wherein the signal added to the original color image is a 
periodic signal which is greater than extent of the original color image. 

5. A method for operating a machine to add humanly imperceptible information to an original color image data struc- 
5 ture to produce a modified color image; the machine including a processor and a memory device for storing data; 

the data stored in the memory device including instruction data the processor executes to operate the machine; 
the processor being connected to the memory device for accessing and executing the instruction data stored 
therein; the method comprising: 

operating the processor to obtain an original color image data structure defining an original color image in- 
cluding a plurality of original image locations each indicating an original color; 
operating the processor to obtain periodic signals indicating information: 

operating the processor to add the periodic signals indicating information to the original color image by mod- 
ifying a plurality of the original image colors in original image locations to produce an interim color image data 
structure; the interim color image data structure being a version of at least a portion of the original color image 
data structure and including a plurality of interim image locations each indicating an interim color; an original 
image location in the original color image data structure and an interim image location in the interim color data 
structure indicating the same image location in each image and being relerred to collectively as respectively 
paired image locations; each one of the respectively paired image locations indicating an image color; 
operating the processor to perform a comparison operation to produce a quantitative perceptual difference 
measurement at each of the respectively paired image locations: the comparison operation using a model of 
human perception to quantitatively represent the perceptual differences as perceived by a human viewer of 
the original and interim color images; 

operating the processor to attenuate the periodic signal to produce a non-zero modified periodic signal when 
the difference measurements between the respectively paired image locations exceed a threshold indicating 
an acceptable perceptual difference between the original and interim color images; and 

operating the processor to perform an iteration operation iterating the adding, comparison and attenuating 
operations using the modified signal until the perceptual difference measurements between the respectively 
paired image locations do not exceed the threshold; the iteration operation producing a modified color image 
data structure including the humanly imperceptible information and indicating a perceptually acceptable signal 
difference between the original and modified color images. 

6. A method for operating a machine to automatically transform geometric properties of an acquired image version 
of an original image to match geometric properties of the original image; the machine including a signal source 

35 connected for receiving signals, a processor and a memory device for storing data; the data stored in the memory 

device including instruction data the processor executes to operate the machine; the processor being connected 
to the memory device for accessing and executing the instruction data stored therein; the method comprising: 

operating the processor to obtain from the signal source an acquired image data structure defining an acquired 
image; the acquired image being a version of an original image and additionally having embedded signals 
therein not included in the original image; the embedded signals having predetermined geometric relationships 
with respect to each other; 

operating the processor to obtain geometric constraint data indicating expected geometric relationships about 
the embedded signals in the acquired image; 

operating the processor to determine the predetermined geometric relationships of the embedded signals in 
the acquired image; 

operating the processor to compute geometric differences between the acquired image and the original image 
using the predetermined geometric relationships of the embedded signals in the acquired image and using 
the geometric constraint data indicating the expected geometric relationships of the embedded signals, 
operating the processor to transform geometric properties of an acquired image version of an original image 
to match geometric properties of the original image using the geometric differences computed between the 
acquired and original images. 

7. A method according to claim 6, wherein the embedded periodic signals include at least two embedded, preferably 
55 sinusoidal, signals. 

8. A method according to claim 7, wherein determining the geometric relationships of the embedded signals in the 
acquired image using the geometric constraint datas includes 
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operating the processor to produce an amplitude spatial frequency spectrum image of the acquired image; 
operating the processor to determine locations in the amplitude spatial frequency spectrum image of local 
peak power concentrations representing the embedded periodic signals in the acquired image; and 
operating the processor to compute the geometric relationships among the embedded period signals using 
the geometric constraint data and the locations of the local peak power concentrations. 

A method for operating a processor-controlled machine to decode periodic signals embedded in an image; the 
machine including a processor and a memory device for storing data: the data stored in the memory device including 
instruction data the processor executes to operate the machine; the processor being connected to the memory 
device for accessing and executing the instruction data stored therein; the method comprising: 

operating the processor to obtain an image data structure defining an image; the image including a set of n 
periodic signals embedded therein; the set of n periodic signals having geometric relationships with respect 
to each other; 

operating the processor to obtain geometric constraint data indicating expected geometric relationships about 
the set of n periodic signals embedded in the image; 

operating the processor to identify a plurality of local peak power concentrations in a two-dimensional spatial 
frequency power spectrum of the image; 

operating the processor to eliminate a subset of the local peak power concentrations that do not satisfy the 
geometric constraint data indicating expected geometric relationships about the set of n periodic signals em- 
bedded in the image; a remaining subset of the local peak power concentrations being determined to be 
candidate local peak power concentrations indicating spatial frequencies and directions of the set of n periodic 
signals; and 

operating the processor to compute a best linear mapping between sets of n candidate local peak power 
concentrations in the image and a set of n periodic signals indicated by the geometric constraint data; the best 
linear mapping identifying the spatial frequencies and directions of the set of n periodic signals embedded in 
the image. 

A method according to claim 9, wherein the image is a color image; wherein the set of n periodic signals embedded 
in the color image are embedded in at least one of a plurality of opponent color image representations of the color 
image; and wherein identifying the plurality of local peak power concentrations in the image includes coverting the 
color image into the plurality of opponent color image representations and identifying the plurality of local peak 
power concentrations in the at least one of the plurality of opponent color image representations. 
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(54) Method for embedding signals in a color image 



(57) A signal component (e.g., a period signal com- 
ponent such as a sinusoid) is added (210) to an original 
color image to form a modified color image, with the goal 
of embedding the signal in the image so that it is imper- 
ceptible to a human viewer. A comparing operation 
(240) uses a model of human perception to measure the 
perceptual difference between the original and modified 
images, identifying local areas of the modified image 
where the signal difference exceeds a threshold, indi- 
cating that the signal is perceptible to a human viewer. 
Using the perceptual difference measurement data; the 
signal is attenuated in the identified local areas that in- 
dicate a perceptually unacceptable difference, and this 
modified signal component is then added (210) to the 
original color image in a next iteration. Perceptual dif- 
ference measurement and signal attenuation are iterat- 
ed until the comparison operation determines that the 
signal difference is perceptually acceptable. The tech* 
nique takes advantage of the fact that, although the spa- 
tial frequencies of the embedded signals are well within 
the range of spatial frequencies to which humans are 
normally quite sensitive in the luminance (black-white) 
vision channel, this sensitivity does not extend to the 
color vision bands. In an illustrated embodiment, a set 
of sinusoidal signals that form a grid are added to the 
color image; location (i.e., decoding) of the sinusoids, 



which does not require the original color image, allows 
computing a geometric mapping from an image with the 
embedded signals to the original image. 
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